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Abstract 

We revisit a two-component inflaton model with a turning trajectory in the field space, where 
the field slowly rolls down along the trajectory. We consider the case when the effective mass in the 
direction perpendicular to the trajectory, namely the isocurvature direction, is either of the same 
order as or much larger than the Hubble parameter. Assuming that the turning angular velocity 
Q_|. is small, we compute analytically the corrections to the power spectrum of curvature perturbation 

caused by the mediation of the heavy isocurvature perturbation, and compare our analytic results 
with the numerical ones. Especially, when M 2 S S> H 2 , we find that it is proportional to M~^. 
' This result is consistent with the one obtained previously by an effective field theory approach. 
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I. INTRODUCTION 



Inflation [j], |2[ has been recognized as the most competitive model for the early universe, 
which can not only solve basic problems of big-bang cosmology, say, the horizon and flatness 
problems and overproduction of topological defects, but also give a natural explanation for 
the primordial fluctuations that account for the observed CMB anisotropy and large scale 
structure of the universe. 

The simplest model of inflation is driven by a single scalar field with a canonical kinetic 
term and sufficiently flat potential, which generates an almost scale-invariant, highly Gaus- 
sian curvature perturbation [2(J . But one can still consider a single-field model with complex 
kinetic term like k- inflation [18J or DBI inflation [Toj. The WMAP has accurately confirmed 
the existence of primordial curvature perturbations with a nearly scale-invariant spectrum 
of order 10~ 9 , but at the same time indicated the possible existence of non-Gaussianity that 
cannot be explained by a single- field slow-roll inflation with canonical kinetic term [3j. 

Multi-field inflation is another natural generalization of the simplest single-field inflation. 
It is motivated partly because of some theoretical considerations, like string landscape Q, 
and partly because it can lead to more abundant phenomena. For example, N-flation or 
assisted inflation [9[ shows that N scalar fields move collectively to sustain inflation to an 
e-folding number required, even if each individual field is unable to drive inflation. Also one 
can consider multi-field k-inflation or DBI inflation 

As a case study, two-field inflation can bring us fruitful new properties while maintaining 
geometrical intuitions in the field space. In slow-roll paradigm, we can classify two-field 
models into different categories by their masses, velocity in field space, etc.. The simplest 
case is when both fields are almost massless. In this case the effect of the curvature of 
the trajectory in field space can be neglected, implying that the equations of motion can 
be effectively decomposed into the one along the trajectory and the other orthogonal to 
the trajectory. One then call the field along the trajectory the "adiabatic" component 
and the one orthogonal to it the "entropy" or "isocurvature" component. The curvature 
perturbation (literally speaking the curvature perturbation on comoving or uniform density 
slices) is directly related to the adiabatic component of the field perturbation (on flat slices), 
and is affected by the entropy perturbation via "transfer function" , which can be calculated 



in super- horizon e ra Il4-|l7l| . Or, one can calculate directly the power spectrum by using 



the 5N formalism 123-26 



Another important case is when the trajectory is along a deep valley, that is, the effective 
mass along the trajectory is sufficiently small so that the field is in slow-roll motion, while 
the effective isocurvature mass perpendicular to the trajectory is heavy, M c fj ^> H . In such 
a model, the behavior is more or less like that in a single-field theory since the isocurvature 
perturbation quickly decays out. Namely, integrating out the heavy field, the system can 
be described by an effective single field theory. In this effective field approach, it has been 
recently claimed that the effect of the heavy field can be absorbed into a correction in the 
sound speed which is inversely proportional to the squared effective mass of the isocurvature 



component |36H38j . However, because of some non-trivial assumptions in this approach it 
is desirable to compute the correction without resorting to the effective field approach, and 
justify or falsify it. 

The case that adjoins these models is called "quasi-single field" inflation by [341 ]. in which 
the effective mass of isocurvature mode M e g is of the same order as the Hubble parameter 
H . Assuming for simplicity that the field is approximately in a circular motion with a 
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small and constant angular velocity, one can use the in-in formulism [5H7J to calculate the 
2-point function of the curvature perturbation by quantizing the curvature and isocurvature 



perturbations as free fields, with the coupling terms as interaction vertices [19|, |35[. They 
calculated the correction to power spectrum in M eS ~ H case, and found it proportional 
to the angular velocity squared with a M e g- dependent coefficient of order 1, which was 
estimated by numerical methods j3o| . 

The aim of the current paper is to apply the method developed in the "constant turn" 
quasi-single field inflation in |34j to calculate analytically the coefficient of the corrections 
caused by the mediation of isocurvaton to the power spectrum. Especially for the case when 
the effective mass of the isocurvature mode is large, M 2 S ^> H 2 , we examine the result 



obtained by the effective field theory approach [36|, |38|, which is found to be consistent with 
our analytic result. 

This paper is organized as follows. In Section [TT] we briefly review the quasi-single field 
inflation and the in-in formulism used to calculate the power spectrum of curvature pertur- 
bation, and calculate its coefficient analytically . In Section III II we extend our calculation 
to the case of a heavy isocurvature mode, and compare our result with the one obtained 
by a totally different method, namely, the effective field approach in which the heavy field 
is integrated out from the beginning. We conclude our paper in Section HVl Some detailed 
calculations are spelled out in Appendix. 



II. LIGHT ISOCURVATURE MODES: QUASI-SINGLE FIELD INFLATION 

In this section we give a brief review of a simple two-field inflation model: the quasi-single 
field inflation. Using the in-in formalism, we derive corrections to the power spectrum of 
the curvature perturbation due to interactions with the isocurvature mode. 

Quasi-single field inflation describes a segment of time during generic inflation when the 
curvature field is massless and undergoes a slow-roll trajectory, while the isocurvature mode 
has mass of order H. To be specific, we consider a motion along an arc with radius R. 
Then we can naturally decompose the field into the curvature mode along the tangent of the 
circle, R9, and the isocurvature mode given by the radial field a. We call a the isocurvaton. 
The Lagrangian of such a system is 




~(R + a) 2 g^d,9d u 9 - l -g^d,ad„a - V st {6) - V(a) 



(1) 



where V ST (9) is a usual slow- roll potential, V(a) is a potential that forms a circular valley 
and traps the isocurvaton at a = cr . 

The equations of motion for classical trajectory is then 

3M 2 H 2 = l -R 2 9 2 Q + V + V SI , (2) 

-2M 2 H = R 2 9 2 , (3) 
= R% + 3R 2 H9 + V£, (4) 
= a + 3H& + V - R 2 9l , (5) 

where R = R + <Jq. Now &q = gives us V ; (ao) = R9 2 . This relation shows that the field is 
not exactly at the bottom of the valley of the potential V, but is slightly shifted from it to 
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provide the centripetal force of circular motion. In what follows we omit the argument o"o of 
V, V", V'", etc., unless there is a chance of confusion. We denote the "rotation speed" by 9 
as in |15| with a minus sign because of the different orientation. In the effective-field-theory 
4l| . they used r]± = 9/H instead. We define the slow- roll parameters along 
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We will not define the "slow-roll parameters" perpendicular to the trajectory since the field 
does not slowly roll in this direction, and since they can be represent by 9 and V". 
Now let us consider the perturbation to 9 and a in the spatially flat gauge where 



a 2 {t)8, 



and 



9(t, x) = 9 (t) + S9(t, x), a(t, x) = a + 5a(t, x). 



(7) 



(8) 



In this gauge in the leading order in the slow-roll approximation, the gravitational effect can 
be neglected, and we obtain the Hamiltonian density, 

(9) 
(10) 

(11) 
(12) 

where and below we omit the subscript from the background quantities. It is known that 
the conserved curvature perturbation on comoving slices TZ^ is given in terms of the field 



Ho — 


a 3 \-R 2 59j 2 
[2 




+ lsa 2 + ^- 2 (dM) 2 + 1 -M 2 s 5a 2 




—C2a 3 5aS9, 


c 2 = 2R9\ 




n = 


-a 3 R5a59 2 


- a 3 9595a 2 + 


aR5a (d i 59) 2 + —V"'5a 3 , 
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M* G = 


V" + 39 2 , 





fluctuation in the flat slices along the trajectory 59 as [2l|, 122 

H59 

n r = 
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(13) 



For the "constant turn" case, c 2 and M 2 S are both constants. M e g is the effective mass of 
the isocurvature field. We split the Hamiltonian into the free part T-Lo and interacting parts 
^2, Next we will treat I-L2 as interacting vertex to the free Hamiltonian "Ho to calculate 
the power spectrum, while this vertex is proportional to 9q/H and depicted in Fig. [TJ ^3 
is only useful when turning to investigate the bispectrum which is not investigated in the 
present paper. This method is valid only when the interacting Hamiltonian is small, i.e., 
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< 1 



(14) 
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FIG. 1: The second order interacting 
vertex %2 which is proportional to 0/H. 



FIG. 2: The leading order to the 2-point 
function of curvature perturbation, me- 
diated by a isocurvature field. 



Then the perturbation theory can be used, and there is no constraint on the value of M c g ~ 
V": It maybe small, of order H, or large. The terminology "quasi-single field" denotes the 
case M e fj ~ H which we review in this section. The case when M e g ^> H will be studied in 
Section [TTT] by the same approach. 

In the interaction picture, we quantize the Fourier components 58^ and 5a^ of the free 
fields 58 1 and 5a 1 , 

58^ = MkOk + wlk a tk> (15) 
K = v k b k + v\bl k , (16) 

where and 6k are the annihilation operators of 5a and 58, respectively, and each one 
satisfies the canonical commutation relation, 

k, = (2n) 3 5 3 (k + k'), [6 k , = (27r) 3 <5 3 (k + k'), (17) 

with all the other commutators equal to zero. The mode functions, «k and v^, satisfy the 
linear equations of motion followed from the Hamiltonian T-L , 

< - -< + k 2 u k = 0, (18) 

T 

2 M 2 
v » _ - v > + p Vk + —f-v k = 0, (19) 

where r is the conformal time defined by dr = dt/a(t), the prime denotes the derivative 
with respect to r, and we have approximated the background by a de Sitter spacetime. 

The solutions to (1T81) and (fl"9l) are given by linear combinations of the Hankel functions of 
first and second kind. Requiring that the solutions approach those in the Minkowski positive 
frequency functions, or those for the Bunch-Davies vacuum [§[, in the limit k ^> aH, 

TT 

i?w k , v k i^=Te~ ik \ (20) 
V2k 

we obtain 



u k = — ^(l + ikr)e- lk \ (21) 



RV2k 3 



and 



Vk = - ze ^+m^H(-rf 2 Hl 1 \-kr), for M 2 S /H 2 < 9/4, (22) 
where i/ = y/9/4 - M^/H 2 , or 



^ k = -ze-i^f^iJ(-r) 3 / 2 ^ ) (-^), for M 2 S /H 2 > 9/4, (23) 
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where fi = ^/M 2 S /H 2 - 9/4. 

Next we consider the interaction terms as perturbations to "Ho- To calculate the power 
spectrum up to tree-level diagrams, we just need to consider the l-i 2 term which describes 
the correlation between the curvature and isocurvature modes, 



Hi = J d 3 x H{ = -c 2 a 3 J • 
The two-point function of 59 2 can be calculated by the in-in formulism as 



(24) 



(59 2 ) = (0| 



Texp (i / dt'H^t' 



to 



Texp 



dt'H^t'] 



to 



|0), (25) 



(0\T[l + i / dt'H^t')- / / dhH^H^t 



to 



to 



to 



■ T 1 + i 



f dt'H^t') - j dh f dttHftjHrfo) 

to J to J to 



59j (t) 
|0>, 



<0|<W?|0> 



'to 

- 2 Re 
+ ■•• . 



dh / dhiOlHj^) 59]H I {t 2 )\Q) 



Uto 



dti / dt 2 (0\59j iyj(ti)i?/(i 2 )|0) 



to 



(26) 

(27) 
(28) 



To calculate the terms sandwiched with the free field vacuum, we use the normal ordering. 
After that only fully contracted terms survive. The correction from a mediating isocurvature 
perturbation is depicted in Fig. [2j This gives the correction to the power spectrum of 59 as 



(27r) 3 5 3 ( Pl + p 2 ; 



±C{y) 
i? 4 p\ 



3 ' 



(29) 



with a factor C(v) which only depends on the effective mass M e g and rotation speed 9 2 via 



9 M e 2 ff 



4 H 2 



M 2 S = V a <j + 39 2 



(30) 



CM 



7T 



dxx- 1/2 H^(x)e ix 



POO POO 

2Re / dx 1 x^ 1/2 Hl 1 \x 1 )e- ixi dx 2 x 2 1 ' 2 (H^\x 2 ))* • 

JO Jx 1 



(31) 



Applying the above to Eq. ( TT3]) . the power spectrum of the curvature perturbation on co- 
moving slices is given by 



V, 



H 4 



n 



4it 2 R 2 9 2 



1 + 8C(i/) ( -| 



(32) 
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The calculation of C{v) for real v has been done in 34j by numerical integral. It is our main 
task to calculate C(v) in the case when v is pure imaginary in Section ITO1 

As a warmup, we do the calculation of C{v) for < v < 3/2 analytically in this section. 
In this case, we first write 



C{y) 



71 



(A - B), 



(33) 



where A and B are the first and second terms, respectively, of Eq. (1311) . and calculate the 
integrals A and B piece by piece: We use the asymptotic from of the Hankel function at 
i C 1 for the integration from to 1, and the asymptotic from at x > 1 for the integration 
from 1 to oo. We can see that, the main contribution to the integral is from < x < 1, i.e., 
the infrared era when the wavelength of the mode has been stretched out of horizon. 

The asymptotic behavior of the Hankel function in the ultraviolet limit is in the form 
of an ordinary wave, but that in the infrared limit needs some caution since there is a 
spurious discontinuity and divergence at integer v. To avoid it one should take more terms 
of asymptotic expansion into account instead of keeping only the leading order term. We 
defer the detailed discussion in Appendix |A] Here we refer to the result given by Eqs. (lAlj) 
and ( T£2j) , 



A 



dx x- 1/2 H^\x)e l 



E 

k=0 



l + i COt V7I 



2 2k k\ \2 v Y{k + u + 1) J 



dx aT 1/2 #W(xy 



dx x^-VV* 



i2 v 



sin VTiY{k — v + 1) J 



1 ^ poo [W i(2x-i/7r/2-7r/4) 

dxx 2k ^ 2 e lx \+ / dx 



x 



E 

k=0 



(~) fc / 1+ZCOt^TT 2k+u+l/i 



2 2k k\ [2»r{k + v + r 



7 (2A; + v + 1/2, -i) 



sin i>7iT{k — v + 1) 



^2k -v+ 1/2, -i) \ + J - e - i 5(^^)(7r + zEi(2i)) 



(34) 



where some of the integrals are expressed in terms of the incomplete gamma function, 

7 ( s ,p) = f P t^e^dt. 
Jo 



(35) 



The incomplete gamma function is defined originally on the real s axis, but can be ana- 
lytically continued to the complex plane with a branch cut from — oo to 0. For a given 
wavelength, the third term in (134|) is the integral from the ultraviolet era x > 1, which 
contributes little to the final result. 

The evaluation of B is a little tricky since there is an overall time ordering operator. 
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Keeping this in mind, we split the integrals on two clXGS cLS 



B = 2Re I dx x xf /2 H^( Xl )e- ix [ dx 2 x~ 1/2 H { u 

JO J xi 

+ 2ReJ dx l x^ ,2 H^\x l )e- lx J dx 2 x^ l,2 H^\ 



2) (x 2 )e- iX2 



+ 2Re 



dx\ x 1 H^\x\)e~ 



POO 

/ dx 2 X2~ l/2 H^\x 2 )e 



(36) 

(37) 
(38) 



As before, the second line (137|) can be neglected since it is small compared to the first line 
( |36|) . while the third line ( )38|) can be discarded directly if one adds a small imaginary part 
ie to x and do the integral in the ultraviolet limit. Now the double integral in the first line 
( 13"6"|) is equivalent to the integration over the upper triangle bounded by the x 2 -axis, x 2 = 1 
and Xi = x 2 . Since the integrand is in the form 2f(xi)f(x 2 ), it is in particular symmetric 
with respect to the interchange of X\ and x 2 . Hence the integral is equal to a half of that 
over the square bounded by X\ = x 2 = and x\ = x 2 = 1. The result is then given by the 
square of the integral /„ dxf(x), that is, 



B 



Re 
Re 



[ dxx 1/2 Y u \ 
Jo 



x)e 



E 

_k=0 



cot V1X 



2 2k k\ \2T(k 
2" 



Re 



sin i/7iT(k — v + 1 

(_)* /(_^2fc+v+l/2 cotz/7r 



E 



2 2k k\ V 2 v Y(k 



dx x ™+»-V*e- lx 
dx i/2 e — N 

T(2k + u + 1/2,0, i) 



k=0 

sin v-KY(k — v + 1 



T(2A;-i/ + l/2,0,f 



-i 2 



(39) 



In summary, we obtain an approximate analytical expression for C{y) by combining (|33|) . 
( )34|) and ( )39|) . It is depicted in Fig. [3] in which the sum over k is truncated at iV = 10 3 . 
We see that our approximate formula reproduces the one obtained numerically in 34J fairly 
well. In fact, as discussed in detail in Appendix [A] the convergence of this series is very fast 
and a truncation at N = 3 is found to give a result which is almost indistinguishable from 
that at iV = 10 3 . 



III. HEAVY ISOCURVATURE MODES: INFLATION UNDER EFFECTIVE 
FIELD THEORY 

For an isocurvature mode which has a large effective mass, that is, M c r/H > 3/2, the 
parameter v becomes pure imaginary. Therefore we introduce a real parameter \x by 
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FIG. 3: The blue curve depicts an approximate analytic expression for C{v) given by Eqs. ()34p and 
(|39p . discarding the contributions coming from the ultraviolet part x > 1. The dots are those read 



from the numerical result obtained in 
at k < N = 10 3 . 



34 ]. The summation over k in these equations is truncated 



Then the correction to the power spectrum is written as 



c(aO 



7T 



-(17V 



(a-/3) 



a 



dx x- 1/2 H^\x)e ix 



(3 = 2Re / d Xl x^ /z H^( Xl )e- lXl / dx 2 x~ i/z (#^(x 2 ))V 



(41) 
(42) 
(43) 



Note the existence of an overall exponential suppression factor e _/i7r . This kills almost all 
the contributions in a and (3 except for the one which contains an exponential enhancement 
factor that exactly cancels it. The above can be evaluated numerically as depicted in Fig. HI 
But we can evaluate it analytically in the limit when the isocurvaton is very heavy, i.e. 
|U> 1, which is the case when one expects an effective single- field description holds. 
First we consider a. The integral in a is 



dx x 1/2 H^\x)e 



dx 0C J (X) 



e lx + i I dxx l/2 N itl (x)e lx . 



'0 Jo JO 

The first term on the right hand side can be integrated directly as 



doc oc J [x)e 



r(i/2 + ¥) wfl + VL , 

2<M(_j)i/2-HMr(l + ifj,) 2 1 \4 + 2 ' 4 + 2 ' + ^ 

in/4 



2^ 



(44) 

(45) 
(46) 



In obtaining this we have used the special value of hypergeometric function when the last 
argument is 1, see (IB5j) in Appendix [Bj and relation of gamma function. Next, we can 
calculate the second integral in (l4"4l . There is a possible spurious divergence if we integrate 
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it directly. Let us add a small positive imaginary part z'e to the exponent to make it converge, 
i.e. to write the exponent as e lx ~ ex . Therefore, the integral is 



dx x 1//2 A r j A1 (x 



cot z/x7iT(l/2 + ifi) 
2^r(l + i/i)(-2ze) 1 / 4 +^/ 2 

2^r(l - ifi) 
smz>7iT(l/2 - «'//)( -2z'e) 1 /4-W2 



,1 ml m 



2 7 4 2 

l z/i ; 

4 ~ T ' 4 



-2z'e 



1 z> 1 m . 1 



2z'e 



(47) 



Then this is the situation when the last argument of the Gaussian hypergeometric function 
is large, and (1B6j) can be applied to connect it to 



coth finY 



2r(i + z» 



2*T(1 + z» 



2 i T(--z» 



r - + r - + 



2r(--<Ai) 



z sinh jjirT (1 — z'//) / 1 z'/i 



2^1 



Z/i 1 Z/i 

Y'4 + Y ; 



1 + ifi; — 2z'e 



ifi 1 i\x 
2~'4 ~ T 



; 1 — z/i; — 2ze 



Now the divergence disappears, and we can take e — > and use (IBlOj) to get 



2 coth /i7iT 



Z/i 



Z/i 



2 l T 





-t 


HI 













z sinh fiirT 



1 Z/i 

4 ~ y 



3 Z/i 

4 ~ ~2 



Next we take fIBip to cancel these gamma functions. We are left with 



— COth /17T + 



sinh /in 



(48) 



Therefore we have that 



a 



7T 



3 /i7r/2 



V2 



sinh /i7r 



+ z 



V^coth 



/i7T 



(49) 



This is an exact result, with no approximations. We see that when the effective isocurvature 
mass is large, this goes as 

2 

a — > — coth 2 /i7r ~ 1. (50) 

7T 



Together with the factor e _M7r in (141]) . we now that the contribution from a is exponentially 
suppressed. 

Now we move to the evaluation of /3, 



(3 = 2Re / dx 1 x; 1/z H^'{x 1 )e- lX1 / dx 2 x~ 1,2 (H^ {x 2 ))* e~ ix \ 



(51) 
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Take the asymptotic behavior ( 1C5|) we have got in Appendix [Cj we can see 



4e^ 

71 /i 



Re 



exp 



4/i 



exp 



.£l VK/i 



4/i 



2£ 2 



(52) 



which contains a double integral. We first deal with the integral of x 2 . As we have done in 
calculating a, we can take the Taylor series of the exponent for x 2 and the integral of X2 can 
be done term by term to get 



j -yt-1/2 
uX 2 X 2 



exp 



_£2 e i*/4_ ia . a 



4/i 



s^r^" i/2 (i eW4+ -)' 



(53) 
(54) 



fc=0 

00 

£ 

fc=o jfe! I fc + 



2/i 



X < X 



k+l/2-in 
1 



4/i 

1/4 



Xi 



fc+l/2— tu 
- 2/1 



4/i 

'_U/4 



,] 2 F 1 (l,^ + 2A;-z>;^ + A;-i/i; 



\3/4 



^1 



4/i 



4/i 



[(- 



- 4/i]^ ( i, I + 2k - i m 3 - + k- ip t 



(55) 



Again, integrating x 2 from to 00 in (154]) will give us a summation of infinities, which is 
unphysical since we know that in (j53p there is an exponential suppression which will converge 
the integral. Thus, we take a UV cutoff at y\, and set y\ — > 00 after all the integrations and 
summations are done. Thus this integral can be expressed by its lower bound X\ and upper 
bound yi as above via the hypergeometric function 2 F\. We make use of the asymptotic 
behavior of this function when /i ^> 1 in Appendix \B\ From (IB9j) . we see that when \x is 
large and k <C /i, the hypergeometric functions with different fc's are the same, and we can 
resum the polynomials to get back to an exponent, which yields 



1 



1/2-iu 
Vi exp 



, 1/4 



V\ U + 



2/1 



4/i 



.r 



1/2-i/i 
1 



exp 



xl/4 



■ii h + 



4/i 



(56) 



Now we can take the limit of integration upper bound ?/i to be infinite. This kills the first 
term in fl56l) . Then, substituting ( 156]) into fl52|) . 



4 /- 00 
/3 = — Re % I dx\ exp 

This integral can be expressed by the error function, 

2 

7^ 



V2xj 
4/i 



2ix\ 



erf (z) 



'dt 



(57) 



11 



C(n) 




FIG. 4: The red curve shows C(/i) given by our analytical approximation (|6U|) . which is also the 
one derived from the effective single- field approach (|62p . The black dots are the result of numerical 
integrations of (|41[) . There is a small deviation when \x ~ O(10) 



as 



f3 = 2 7 / 4 ^=e- 2 ^Re[ierf(z2 3 / 4 /i)]- 

A/TT/i 3 



(58) 



We can see that it depends on a power-law series of \x when \x is large. To be specific, let us 
can take the Taylor series of j3 around infinity to get 



= -2- 



7T/r 



(2Jfe 



fe=i 



2 5 */V 



(59) 



This is the only term that contributes to the final correction to the power spectrum, and we 
see that the leading term is proportional to /i~ 2 (thus to M 2 S ) when the effective mass of the 
isocurvature mode is large. We take this leading term and go back from the combination of 
(152"]) . (14"T]) and (1591) to obtain our final result 



7? 



P. 



(0) 



1 + 2 



H 2 
MS 




(60) 
(61) 



As before ~ denotes neglecting the terms of higher power. Note the dependence on \i. 
It is not an exponential suppression, but a power-law one. We reached (1611) under the 
approximation that \x is large as well as 9/H is small. As depicted in Fig. HI the analytic 
approximation ( l60i) reproduces a numerical evaluation of the formula ( 14T1) well. This is our 
main result. 

To compare this with the one obtained by the effective-single-field approach, let us refer 



to the result in [36|, |38|, |40|, |4l| where the two-field case with large isocurvature mass is 



considered. After integrating out the heavy field, the power spectrum is corrected in terms 
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of the effective sound speed as 




°. =l+v*T IF > (62) 



where to leading order M 2 S = V aa — 6 in our notation. We first note that the definition of 
effective mass M 2 S is different from ours by 48, which is not important since it is assumed 
to be small when the perturbative in-in formulism can be applied. Besides, we also see that 
our result ( 1601 . which gives the power spectrum a correction factor corresponding to c~ l , is 
exactly the same as ( )62|) at leading order when 9 2 is small. So we have shown the equivalence 
of both methods applying on the same model. 



IV. CONCLUSIONS AND DISCUSSIONS 



In this paper, we studied a two-field model of inflation in which the curvature mode, 
that is the mode along the inflationary trajectory, is light while the isocurvature mode, 
that is perpendicular to the trajectory, is heavy. To simplify the problem, we focused on a 
situation in which the field is in a circular motion with a small constant angular velocity. 
By computing a process that describes the effect of the heavy isocurvature mode on the 
curvature perturbation 2-point function by the in-in formalism, we derived a correction to 
the curvature perturbation power spectrum due to the heavy mode intermediation. The 
result tells us that the correction is proportional to (8/M c g) 2 , which is exactly the same as 
the one obtained by the effective single-field approach by integrating out the heavy field. 

Here let us reconsider our assumptions, and discuss the validity and limitaions of both 
methods. We already know that only if 



H 2 



< 1 



(63) 



one can use perturbation theory in the in-in formalism which we employed. On the other 
hand, in the case of a sharp turn with large angular velocity, even for a very short time, 



39, 45 



that such cannot be described by 

The other assumption is M 2 S ^> H 2 which we use to 



our method fails. It is demonstrated in 
an effective single-field action either, 
get our final analytic expressions. Although this is not necessary for the use of the in-in 
formulism, it enabled us to obtain a simple analytic result and to easily compare this with 
the effective single-field approach in which the heavy modes are integrated out. 

Some papers discussed the requirements of effective single field theory to hold besides the 
isocurvature mode is heavy. The main constraint is called the adiabatic condition 



40 



dt 



log 9 



cff • 



or 



dt 



iog(c; 2 - 1) 



cff • 



(64) 



In our constant turn case these conditions are satisfied since we have 6 = 0. It is shown that 
the decoupling of the heavy field by integrating it out is valid only when slow-roll condition 
is preserved (42|] . In [39| the authors claim that to keep slow- roll in the classical trajectory 
one should impose 

V" 9 

H»H- < 65 > 
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FIG. 5: This is the main contribu- 
tion to large non-Gaussianity eval- 
uated in [34f | . 



FIG. 6: These are other possible subleading contributions to the 
non-Gaussianity in quasi-single field inflation. 



which is satisfied for \i ^> 1. On the other hand, for v — > 3/2 which corresponds to the 
case when not only the inflaton but also the isocurvaton have their masses negligible, there 
does not exist even an appropriate classical quasi- single-field trajectory. This is probably 
the reason why we see the divergent behavior in Fig. [3] at v — > 3/2. In this case, one 
should abandon the effective single-field description, and turn to other methods used in the 
standard slow-roll multi-field case as mentioned in Section [U 

It is important to study the implications of quasi-single field models in realistic situations. 
It surely represents a kind of turning process, but it must be adiabatic. If the trajectory is 
straight before and after the adiabatic turning, the power spectrum (1611) is only modified 
during the turning process, which implies a feature at a characteristic scale determined by 
the time-dependence of 6. It will be interesting to study this "adiabatic turn" quasi-single 
field inflation in more detail. For instance, we know this scale-dependence can only be 
detected if the turning occurs not too late at the inflationary stage. In hybrid inflation (27[ , 
inflation ends by a sharp turn in the field space and the mode that was isocurvature during 



inflation turns into the waterfall field and becomes important [28H32| . In particular, recently 



it was pointed out that a spiky feature may appear in the spectrum and bispectrum for a 
certain parameter range of hybrid inflation 33j. The waterfall process is not adiabatic. Thus 
our method cannot be applied. Nevertheless, it is interesting to see if an extension of the 
present method based on the in-in formalism can shed more light on these cases of hybrid 
inflation. 

Finally, let us comment on possible developments from our current work. We notice that 
our result flBTT) is exact only at leading order in the limit fi — > oo. Beyond this limit we 
can calculate higher order correction in 1/M e g to compensate possible deviations from the 
leading order result. An interesting possibility emerges when we consider next-to- leading 
order terms. We know in the EFT approach, corrections to the power spectrum from the 
isocurvaton come in only via V aa , which means the next-to-leading order correction should 
be of order 1/Vj 2 ^ ~ 1/M* S . On the other hand, in our approach, the effect comes from 

= ^M^/H 2 — 9/4 ~ M cS , thus naively the next-to- leading correction will be 1/M^ S . 
Therefore, unless we has a mechanism to cancel all the contributions from the terms with 
odd power indices, we will face this inconsistency. We need further investigations to find 
out the possible implications behind this fact. 

Another access to check this consistency is to consider the higher order correlations, 
specifically the 3-point function of the curvature perturbation. Actually, one of the 
main motivations for developing the quasi-single field approach was to calculate the 
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non-Gaussianity expected in such a model [34|, |43j. It is found that only the vertex given 
by V" , i.e., the last term in (TTTj) corresponding to the diagram shown in Fig. [5] was of most 
importance in contributing to large non-Gaussianities. This is because the self-interaction 
of the isocurvature mode, i.e. V", is possible to be large, while the other vertices in ( ITT]) , as 
are depicted in Fig. El will give a negligible contribution. Under the same assumptions, we 
can also calculate the similar coefficient s{v) and a(v) of non-Gaussian parameter /nl in 
the three-point correlations which also depend on the effective mass, extrapolate it to the 
case with a large isocurvaton mass, and compare it with the non-Gaussianity in effective 



field approach [36|, |38|. It is very interesting to see whether these two results still coincides 
with each other, and if so, whether we can prove the equivalence of the two approaches to 
all orders. 

Note: When this work was nearly completed, we were acquainted by private communi- 
cation that a related work j3 on this issue was also approaching its completion. We thank 



the authors of [44J for synchronizing our submission to the arXiv. 



Acknowledgments 

We are grateful to Yifu Cai for useful discussions on numerical methods. SP thanks 
the hospitality during his visit to Yukawa Institute for Theoretical Physics for the workshop 
"2012 Asia Pacific School/ Workshop on Cosmology and Gravitation" (YITP-W-11-26) when 
this work was initiated. SP also thanks Bin Chen, Qiang Xu, Jiaju Zhang for valuable 
comments on the work. We thank the referee and editor of JCAP, who have pointed out 
the inconsistency between our result in the first version and EFT's, and give us suggestions 
on how to improve the accuracy in our calculation. This work is supported in part by the 
JSPS Grant-in-Aid for Scientific Research (A) No. 21244033, and by the MEXT Grant-in- 
Aid for the global COE program at Kyoto University, "The Next Generation of Physics, 
Spun from Universality and Emergence". SP is supported by the NSFC Grant No. 10975005 
and Scholarship Award for Excellent Doctoral Student granted by Ministry of Education of 
China. 



Appendix A: The Asymptotic Form of Hankel Functions of Real Order 

Here we give several useful expressions for the asymptotic expansion of the Hankel func- 
tions, paying particular attention to the case of integer orders 
We start from the series definition of the Bessel function, 



JJx) 



E 

fc=0 



(-)' 



1 



fc! T(k + u+ 1) V2 



) 



(Al) 



From this definition the Neumann function is defined by 

1 



YJx) = cot vkJJx) 
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COt Z/7T 
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■ (A2) 
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Forz/ > 0, it seems that we can neglect the first term in ( IA2I) since it is small compared 
to the second one when x — > 0. However, both terms suffer from divergence when v is an 
integer, hence the first term cannot be simply discarded. 

To see how the divergence in each term is canceled with each other, let us consider first 
the behavior at v = 0. In this limit, we have 



1 



lim YJx) = — 1 1 + v In ( — 

oo 

fe=l v 

= -J (x)ln^ 

2 1 X t 

— > — m — tor i-^ 0. 

vr 2 



1 - v In - 



z/ln 



1 — vhi 



(A3) 
(A4) 



Thus we obtain a finite result because of the cancelation of 1/V divergences coming from 
both terms. 

Now consider the case when v is an integer. In our case, since v < 3/2, only the case 
v = 1 concerns us. Introducing e = 1 — u, we have 



lim Y u (x) 

v — ^ 1 



£ 

fc=i 



(-) fc /£\ 2 



1 



e7rr(A; + 2) V2 



-r\ 1 



'2 

em — 

.2 



1 



1 (l + eln( 



(A5) 



~% to finite terms. 



enT{k + e y 

= + -Ji(ar)ln(- 

2 

— s> when x — > 0. 

Thus there is no divergence as 1/(1 — v)- 

Something subtle appears when we want to truncate the summation 
From the calculation above, we see that each term in the series is divergent, although the 
residue is suppressed exponentially by l/r(A; + l)T(k + u+ 1) or l/T(k + l)T(k — u+l) when 
k is large. Expanding (IA2j) . we find that for an integer u, the divergence of the k + -th term 
of positive series, x 2k++I/ , is canceled by the divergence of the k~-th term of the negative 
series, x 2k ~ u , where k~ = k + + v. Truncating the k + and k~ series simultaneously at the 
same integer will cancel the divergence at v = 0, but all the other divergences at v > 1 
survive. Similarly, one can also truncated at k + = N and k~ = N + m, where m is a positive 
integer to cancel the divergence at v = m, but divergences at all the other integers remain. 
Thus unless we take the infinite series limit A^ — > oo, all the divergences at all the integers 
cannot be canceled simultaneously. 

Fortunately, however, the residues of these divergences are suppressed by T(N) for large 
N. So in practical contexts this will not cause a problem. For example, if we truncate the 
series at k + = k~ = 1000 to cancel the divergence at v = 0, the one at v — 1 remains. 
However since the residue is suppressed by l/r(1000) 2 , the divergence becomes effectively 
invisible unless we take the value of v exponentially close to v — 1, within the width of 
Au ~ l/r(1000) 2 . This enables us to draw our analytic curve practically smooth as depicted 
in Fig. [3j 
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With the above understanding, we obtain the Hankel functions in the form, 

H^\x) = J v {x) + iY v {x) ->■ iY v (x)\ (A6) 
H {2) (x) = J v {x) - iY v {x) -»■ -iY v {x) (A7) 



where the Neumann function is given by the form (|A2j) . The last arrows above hold only 
for a real v but not for an imaginary u, ^ = z/i for /i > 0. 
As for the ultraviolet behavior at x ^> 1, it is trivial, 
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(A8) 
(A9) 
(A10) 
(All) 



These formulae hold for any complex order v. 



Appendix B: Gamma Function and Hypergeometric Function 

A useful relation which is used in the context is the product of two gamma functions 
whose arguments vary by 1/2: 

r(*)rQ + *) =^r(2*). (bi) 

Next let us refer to the asymptotic expansion of the Gamma function T(z) at z — > oo on 
complex plane, 

1 1 



Viz) ^(z) z ~ l/2 e~ z V2^ \l + ^ + ...}. (B2) 

v ; v ; \ \2z 288^ 2 J v ; 

In this paper the most commonly used argument is like a + i\i where fi ^> 1 and a <C 
Substituting it into ( 1B2|) we have 

T(a ± z/x) « v ^ e ±i[/iln ^ /i+ " (a ^ 1)/4] /i a " 1/2 e^ /2 . (B3) 

Here for simplicity we only preserve the leading term. 

The Gaussian hypergeometric function 2-^1(0, b; c; z) appears when dealing with the inte- 
grals in section IHH 

rw , x T(a + n)T(b + n)T(c) z n 

^ Y{a)Y{b)Y{c + n) n! 

We will use the hypergeometric function in three different limit of z. First, from the defini- 
tion, we have 

,F,(a,ft 7; l)= y r(7 - a ~fl> (B5) 
r(7-a)r(7-/3) 
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Now deal with large z, which can be connected with the hypergeometric function with a 
small argument 1/z as 

r(aW U (<*,/?; 7;,) = r(cW - «) ^ f g _ + g _ p + i\ 



r(7) H7-a) V z 

+ r r ) h-^ ) (--)'^ - 7 + - « + li g -(B7) 

It holds when |z| > 1 and 2 is not a positive real number. 

The other limit is z — > 0. In this paper we are interested in the arguments a — 1, 
6 = 3/2 + 2& — ifi and c = 3/2 + k — ifi. Thus we have 



z n . (B8) 



r - + 2fc-i/x + n r - + A;-2/x 



/ 3 3 \ V 2 7 V 2 7 

2 F 1 h,- + 2k-ifi;- + k-ifi-z\=J2 —7^ \ 7^ \ 

«=°r f- + 2fc-i/xj T l- + k-ifi + n) 

When we take the limit fi — > 00, k <C fi and n <C /x, using ( 1B3I) . we can resum the power 
series to get 

/ o o \ 00 1 

2 Fx M , - + 2A; - i/i; - + A; - ifi; z J « ^ = — . (B9) 

^ ' n=0 

The approximation holds because \z\ < 1 and the main contribution comes from the terms 
with small n («C fi), which makes it possible to push the summation upper bound to infinity. 
When we take z — > 0, we have 

2 F 1 (a,/3; 7 ;0) = 1 (BIO) 
which can also obtained by setting z = in the definition (1B4I) . 



Appendix C: Hankel Function of Imaginary Order 

In this appendix we turn to the Hankel function of imaginary order ifi, focusing on the 
large fi limit, which corresponds to the isocurvaton with a large mass. For an imaginary 
order, in contrast to the case discussed in Appendix El we cannot neglect J^. We have 



H w ( ) = V ^ l+tCOtiflTT /X\2fc+»M 

n * W ■ " 2^ fc! rOr+l + zV) V27 

' <£ \ 2k— ifi 



k=0 



smifin k\ r(k — l + ifi) \2/ 



ti rrt - 1 4-»ii.^ (9) ' ( C1 ) 



We consider the large mass limit of (jCTf . that is, /i 00. First note that in this limit, we 
have 

cot ifin = — i coth fin — > —i, (C2) 



le 



sin i/i7r = i sinh //7T — > — — . (C3) 
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The above expressions imply that the second line of (1C1I) is exponentially suppressed. Be- 
sides, since the large k terms in the summation of the first line of (ICip are also highly 
suppressed by the factor {T{k + l)T(k + 1 — ifj,)) -1 , we expect the main contribution to this 
summation comes from the terms with small k. This permits us to apply k <C fi in the 
Gamma function, and use the approximate form given in (1B3|) . 



F{k + l+in) 



2 7re «(M In „fc+l/2g-/XTr/2 
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Therefore, we can pick out the /c-dependence in the Gamma function, and resum the poly- 
nomial to get 
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